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ABSTRACT 



Metallicity in the Galactic Center: 
The Quintuplet cluster 

Francisco Najarro^, Don F. Figer^, D. John Hillier^, T. R. Geballe^, Rolf P. Kudritzki^ 

naj arroOdamir . iem .csic.es 
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Q^i We present a measurement of metallicity in the Galactic center Quintuplet 

Cluster made using quantitative spectral analysis of two Luminous Blue Variables 

(LBVs). The analysis employs line-blanketed NLTE wind/atmosphere models fit 

to high-resolution near-infrared spectra containing lines of H, He I, Si ll, Mg ll, and 

Fell. We are able to break the H/He ratio vs. mass- loss rate degeneracy found 

""q^iI in other LBVs and to obtain robust estimates of the He content of both objects. 

Q i Our results indicate solar iron abundance and roughly twice solar abundance 

^ I in the a-elements. These results are discussed within the framework of recent 

d [ measurements of oxygen and carbon composition in the nearby Arches Cluster 

and iron abundances in red giants and supergiants within the central 30 pc of 

, the Galaxy. The relatively large enrichment of a-elements with respect to iron 

lO ! is consistent with a history of more nucleosynthesis in high mass stars than the 

00 

I I Galactic disk. 

rn 



Subject headings: Galaxy: abundances - stars: abundances - stars: individual 



O I (Pistor Star, FMM362) - infrared: stars - Galaxy: center 

O 

^ . 1. Introduction 

Elements heavier than hydrogen and helium ( "metals" ) are primarily created by nucle- 
osynthesis in stars. Metals are important ingredients in many astrophysical processes such 
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as radiative cooling, and mass-loss during star formation and at all stages of stellar evolu- 
tion. They also play a fundamental role in stellar evolution through their influence on stellar 
opacities, and represent a historical record of galactic chemical enrichment via stellar winds 
and supernovae ejecta. 



seen in stars and gas ( 


Afflerbach et al. 


1997; 


Rudolph et a . 


2006; 


Maciel & Ouireza 1999: 


Fuhrmann 1998: RoUeston et. al. 


2000 


ISmartt et. al. 


2001: 


Luck et al. 


2006 


) . Other galax- 


ies show a similar trend, having highest metal abundances in their nuclei ( 


Urbaneja et al. 



20051 ; iKennicutt et. al.ll2003f ). 



P revio u s wor k ( Frogel et al.lll999 ; Feltzing &: Gilmor^boOO : Carr et aPboOO : Ramirez et al 



19971 . Il999l . I2OOOI ) on the Galactic center (GC) has indicated ro ughly solar stellar metal 



abundances, whereas the analyses of interstellar emission lines ([Shields fc Ferland I Il994 



Maeda et al.ll2000l ) have suggested considerably higher abundances. It is not clear why the 
stellar and gas-phase measurements should differ so greatly. 

The GC contains three dense and massive star clusters that have recently formed in 
the inner 50 pc, the Arche s, Qui ntuplet, and Central clusters. Using quantitative spec- 
tral analysis, iNajarro et al.l ( 20041 ) (Paper I) determined that the WNL stars in the very 
young (2-2.5 Myr) Arches Cluste r have roug : hly solar raetallicities. Being more evolved 



4 Myr), the Quintuple t Cluster (IGlass et al.lll987l . Il990l ; iNaeata et al.lll990l ; lOkuda et al. 



I990I ; iMoneti et al.lll994l ) contains a variety of massive stars, includ ing WN, WC, WN9/0fpe, 



1995; 


Fieer et al. 


1999a 


bf). Two LBVs in it are known, the Pistol Star ( 


Moneti et al. 


1994; 


Cotera et al. 


I1994I: 


Fieer, McLean, & Morris 


.995: Fieer et al. 


1998. 1999c) and FM 


M362 
n rich 


(Fieer et al. 


1999b: 


Geballe, Najarro, & Fieer 


2000), each havine an infrared spectrur 



in metal lines of Fell, Sill, & Mgll. 

In this paper, we use quantitative infrared spectroscopy of the two Quintuplet LBVs to 
make direct determinations of metallicity in those stars. We also use the derived a-elements 
vs. Fe ratio to address the dominance of massive stars on the IMF in this region. 
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Fig. 1. — Model fits {dashed lines) to the observed infrared diagnostic lines {solid lines) of 
the Pistol Star. The forbidden [Fell] line at 1.677 fim was not included in the models. 
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Observational Data 



The data were obtained at UKIRlQ usine; CGS4. The Pistol Star was observed in 
April 1996 (Pa; R~3000), July 1997 (L; R~16000) and April 1998 (H; R~5000 and K 
R~3000). Likewise, spectra for FMM362 were obtained in April (L) and May 1999 (H 
and K), using CGS4 in medium resolution mode (R^5, 000-6, 50 0). The slit width was 0'/6 
for all observations. We used the photometric measurements of iFiger et al.l (119981 ) for the 
Pistol Star, to scale the reduced spectra. For FM M362, given its ph otometric variability, 
we adopted the average value, K=7.30, obtained by lGlass et al.l (119991 ) for the epoch closest 
to our spe ctroscopic observa tions. This value agrees, within t he 0.26 standard deviation 



derive d by lGlass et al.l (119991 ) . with the K=7.50 value adopted by lGeballe. Najarro. fc Figer 



(I2OOOI ) from flux-calibrated spectra. We assume the sam e extinction for both objects and 
adopt the value of Aii-=3.2 derived by iFiger et al.l (119981 ) for the Pistol Star. The reader is 
referred to these papers for a detailed discussion on the reduction of the observed spectra 
and photometry. 



Models 



To model the LBVs and estimate their physical parame ters, we have used CM FGEN, 
the iterative, non-LTE line blanketing method presented by iHillier &: Millerl (119981 ) which 
solves the radiative transfer equation in the co-moving frame and in spherical geometry for 
the expanding atmospheres of early-type stars. The model is prescribed by the stellar radius, 
i?*, the stellar luminosity, L^,, the mass-loss rate, M, the velocity field, v{r) (defined by V^o and 
/?), the volume filling factor characterizing the clumping o f the stellar wind, f(r) (see Sec. l4.2p . 
and elemental abundances. iHillier fc Millerl (Il998l . Il999l ) present a detailed discussion of the 
code. For the present analysis, we have assumed the atmosphere to be composed of H, He, C, 
N, O, Mg, Si, S, Fe and Ni. Given the parameter domain the LBVs are located, the r = 2/3 
radius is located close or above the sound speed, and therefore the assumed hydrostatic 
structure plays no role. Thus, no spectroscopic information abo ut the mass of the object 
can be obtained. The atomic data sources are described in detail IHillier et al.l (|200ll ). Here 
we focus on the model atoms used for our abundance determinations Fell, Mgll and Sill. 
Using the supe rlevel formalisra (NS/ NF, number of superlevels vs. number of levels in the 
full atom, e.g. iHiUier fc MiUeil Il998h we we chose 233/709 (up to 109800 cm'^) for Fell, 
37/50 (up to 119400 cm'^) for Mgll and 35/72 (up to 125000 cm^^) for Sill. The choice 



^The United Kingdom Infrared Telescope (UKIRT) is operated by the Joint Astronomy Centre on behalf 
of the Particle Physics and Astronomy 
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Fig. 2. — Model fits {dashed lines) to the observed infrared diagnostic lines {solid lines) of 
FMM362. 
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of the appropiate packing has been extensively tested in iNajarro I (120011 ) . We will revise the 
importance of this issue for the case of Mg in 14.41 

Observational constraints are pr ovided by th e H, K a nd L-band spectra of the st a rs and 
the dereddened K ma gnitud es from Figer et al.l (119981 ) . iGeballe. Najarro. fc Figer I (120001 ) 
and I Glass et al.l (119991 ) . As in lPaper j . a distance of 8 kpc has been assumed. The v alidity of 
our technique has been demonstrated in iNajarro et al.l (1l999l ) and iNajarro I (l200ll ) by cali- 
brating our method against stars with similar spectral type such as P Cygni and HDE 316285 
for which not only infrared but also optical and UV spectra are available. 



4. Results 



Table 1 gives the derived stellar parameters for both LBVs and Figs. [T] and [2] show 
model fits to the relevant lines in the stars. Theoretical spectra have been convolved with 
the instrumental resolution. We note that given the large number of parameters involved 
in the analysis it is unafordable to perform a full systematic error analysis in the whole 
parameter domain. We rather proceed by estimating the range of values for the main stellar 
parameters which provide acceptable fits to the observed spectra. Once those ranges are 
set, we derive the corresponding abundances and their errors. From Table 1, it can be seen 
that the Pistol Star and FMM362 have very similar properties, with the exceptions of the 
Pistol Star's significantly higher wind density (evidenced by the its stronger spectral lines) 
and its higher He content. The latter may denote a slightly advanced evolutionary stage for 
the Pistol Star (see below). Given the general resemblance of the spectra of the objects, we 
discuss them together. 



4.1. Main Diagnostic Lines and Stellar Properties 

Several spectral diagnostics constrain our estimates of the stellar temperature, and thus 
the ionization structure, in particular the He I (5-4) components near 4.05 fim. If helium is 
predominantly singly ionized, even for the most favorable case with (minimum) cosmic helium 
abundance, the observed ratio of H to He I lines exceeds the expected values by large factors. 
This indicates that Hell must recombine to He I very close to the photosphere, implying 
an upper limit of around 13,000 K for the temperatures of these objects. We find a lower 
limit of 10,000 K for the temperatures, as lower values would require non-detection of the 
He I components. Also the strengths of the He I 1.700 fim and He I 2.112/3 fim lines are very 
sensitive to temperature, so that they appear in emission above 12,500 K and vanish below 
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Table 1. Derived stellar parameters 



Parameter 


Pistol 


FMM362 




1.60 


1.77 




dUb 


350 




1.18 


1.13 


H/He 


1.5 


2.8 


Fe/Feo 


1.1 (0.78) 


1.1 (0.78) 


Mg/Mgo 


2.2 


1.5 


Si/Si© 


1.8 


2.1 


1sA{W-^Mq yr-^) 


2.1 


1.2 




3.0 


1.3 


Ko(kms-i) 


105 


170 


D^o^=log(MKo(R/i?0)'/') 


29.39 


29.38 




22.5 


30.5 


CLi 


0.08 


0.08 


CL2 


2.5 


2.00 


CL3 


2.00 





Note. — p is the exponent des cribing the velocity field, D^om is the modified wind 
momentum (IKudritzki fc Pulsll2000l ) and clumping parameter, the CL are defined in Eq. [H 
M^Edd are the Eddingt on masses. H/He i s the ratio by number, and other abundances are 
relative to solar af t er iGrevesse fc Noeld (119931 ). The Fe abundance relative the value in 



Anders fc Grevessd (119891 ) is given in parentheses (see text). 
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10,500 K. These lower limits on the effective temperature are also consistent with the non- 
detections of the Sin 3s'^3p'^Sy2-3s^^p'^P3/2 2.180 /xm and Ss^Sp^Si/s-SsMp^Pi/s 2.209 /im 
intercombination lines, as they are expected to appear strongly in absorption as soon as 
the temperature drops below 10,000 K. Hence, AT±1500 K are conservative estimates of 
uncertainties for the temperatures of the Pistol Star and FMM362 given in Table 1. 

To estimate the termin al velocities, we make use of th e Fe ll] (serni-forb idden Fell) 
z^Fg/a-c^Fg/a 1.688 fim line dCeballe. Naiarro. fc Fig^boool : iFiger et allEoOsI ) that forms 
in the outer wind and has a weak oscillator strength (gf~10~^). This is because non- negligible 
continuum opacity effects at 4 /im may provide only lower limits if Br a is used. The larger 
Voo derived for FMM362 can be clearly inferred from the width of this Fell] line and the 
obvious overlap at Ba between the Hi and He I components (see Fig. [2]). For the Pistol Star 
(Fig. [T]) the Ba components are fairly well separated. 

Our analysis of wind den sity (M, 0) giv es values of (3 that agree with those inferred i n 
the literature for other LBVs (INajarro Il200ll ) and B-Supergiants (e.g lCrowther et al.ll2006l ). 
The values are fairly well constrained by the shapes of the hydrogen lines, especially those 
of Bo; and B7 which are inconsistent with the same value of (3 for both objects. 

Although the wind density derived for the Pistol Star is much higher than for FM M362, 
the modified wind momenta Dmom = log {MVao\/R/ Rq) i Kudritzki Sz Puld I2OOOI ) of the 
two LBVs are nearly identical (Table 1 ). This result is qualita tively consistent with the 
wind momentum - luminosity relation (IKudritzki fc Pulsl I2OOOI ) which predicts the same 
modified momenta for objects with the same stellar type and luminosity. Note that we 
have used dumping-corrected values of M to compute the modified momenta. If unclumped 
values were assumed, Dmom would be closer to 30.0. Interestingly, the latter agrees very 
well with the averaged modified momentum of AG Car at maximum ( i.e., at simi l ar Tpw 
to our objects) obtained using the values of M, Voo and derived by IStahl et al.l (1200 ll ) 
from fits to Ha. Those authors obtained log(M)~-4.1 for about this temperature, while our 
unclumped values are log(M)=-4.1 and -4.4 for the Pist ol Star and FMM362 r espectively. 
Further, usin g radiation-driven wind models for LBVs, IVink fc de Koted (120021 ) were able 
to predict the IStahl et al.l (1200 ll ) M value for AG Car assuming Voo/Vcsc~l-3 and a current 
stellar mass of SbM^. The same value of K»/Vcsc for the LBVs would imply current stellar 
masses of 27.5Mq for the Pistol Star and AGMq for FMM362. Although these masses should 
be regarded with caution, they are consistent with the Pistol Star being more evolved than 
FMM362 (as inferred from their He/H ratios) and hence having lost more mass during its 
evolution, as indicated by the presence of a nebula around it. 



Compared to the results obtained in iFiger et al.l (Il998l ) by means of non-blanketed 
models, the new blanketed models provide a significant improvement in our knowledge of 
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Fig. 3. — Model spectra showing the effect of clumping parameter CLi (see EqJT]) on diag- 
nostic line profiles of FMM362. Values of CLi are listed in the bottom left panel. 
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the physical properties of these two stars. The degeneracy of the "high" and "low" luminosity 
( Tcfj) solutions for the Pistol Star presented by lFiger et al.l (Il998l ) is broken by the Si ll, Mg ll 
and Fe ll lines, which are clearly more consistent with the "low" solution. We derived for this 
star a luminosity of ~1.6(10^) Lq, an effective temperature of ~11,800 K, and an initial mass 
of 100 Mq. The stellar luminosity is reduced by a factor of two compared with the previous 
estimate, illustrating the importance of the new generation of line-blanketed models. Below, 
we discuss in detail the role of two additional stellar properties that are derived using the 
new models, wind clumping and elemental abundances. 



4.2. Clumping 

Clumping is normally invoked in stellar winds to explain inconsistencies arising between 
p (density) and diagnostics. For a given mass-loss, clumping causes an enhancement of 
processes while leaving unaltered those which depend linearly on p. Further, if mass-loss 
rate and clumping are scaled without changing the M//°-^ ratio, the p-dependent diagnostics 
vary while the recombination lines profiles (oc p^) remain basically unaltered. 

To investigate the clumping we introduce the following clumping law: 



(V-Voo) 



/ = CLi + (1 - CLi)e^ + (CL4 - CLi)e^^ 



where CLi and CL4 are volume filling factors and CL2 and CL3 are velocity terms defining 
locations in the stellar wind where the clumping structure changes. CLi sets the maximum 
degree of clumping reached in the stellar wind (provided CL4>CLi) while CL2 determines 
the velocity of the onset of clumping. CL3 and CL4 control the clumping structure in the 
outer wind. Hence, when the wind velocity approaches V^, so that (V-Vi3o)<CL3, clumping 
starts to migrate from CLi towards CL4. If CL4 is set to unity, the wind will be un clumped 
in the oute rmost region . Such beh avior was already sug gested by lNugis et al.l (119981 ) and was 
utilized by iFiger et al.l (120021 ) andlNaiarro et al.l (120041 ) for the analysis of the WNL stars in 
the Arches Cluster. Recently, iPuls et al. I J2006h also have found similar behavior from Ha 
and radio studies of OB stars with dense winds. Furthermore, o ur clumping paramet r izatio n 
is consistent with results from hydrodynamical calculations by iRunacres &: Owockil (|2002l ). 
From Eq. [Tjwe note that if CL3, and therefore CL4, is n ot considered (CL3 ^0), we recover 
the simpler variation proposed by lHillier &: Millerl (119991 ). To avoid entering free parameters 
heaven we set CL4=1 in all of our investigations, aiming to get an appropriate amount of 
leverage on the amount of non-constant clumping in the outer wind regions. 



CLi: Estimating the wind dumpiness. Figure [3] illustrates the sensitivity to CLi 
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Fig. 4. — Influence of adopted clumping structure (see Eq. [T]) in profiles. Run of clumping 
and line profile sensitivities for different values of CL2 {left panel) and CL3 {right panel), 
compared to observed profiles in the Pistol Star. See discussion in text. 
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of the main diagnostic lines utilized to obtain the degree of clumping in the stellar winds 
of FMM362. For each value of CLi displayed in Fig. [3] the mass-loss rate in the model was 
scaled while keeping M//°'^ constant as described previously. Although there are some lines 
that follow this scaling quite well (Ba, and also B7 and some metal lines not displayed in 
the figure), the He I lines and weak Hi lines react quite sensitively to the absolute degree 
of wind clumping. It can be seen that the He I lines only provide an upper limit to CLi 
(a bit lower than 0.1) and do not react to lower values, but a unique value for CLi can 
be selected by some Hi lines. The Hui4 (Hl(14-6)) line displays the highest sensitivity to 
clumping. Unfortunately, the wavelength interval surrounding this line was not observed 
with sufficiently high S/N in FMM362 and there is a large uncertainty in the continuum 
value, which is critical for estimating CLi. We could make full use of the Hui4 line to 
determine the clumping only for the Pistol Star, where this line is relatively stronger in 
emission (see Fig. [T]). Nevertheless, Fig. Oshows that CLi lies between 0.1 and 0.05. It must 
be stressed that only with a well determined clumping may we address the He/H abundance 
issue (see Sec. 14. 3p . 

CL2 8z CL3: Mapping the clumping structure. The upper panels of Fig. H] illus- 
trate the behavior of the clumping structure for different sets of CL2 and CL3 values, while 
the lower panels display the influence of such behavior on diagnostic lines in the spectrum 
of the Pistol Star. It is evident that for some spectral lines e.g.. He I 2.112 fim, the behavior 
of the profiles with clumping is far from being monotonic. Furthermore, not only do lines of 
different ions react differently to clumping, but also lines within the same ion, e.g.. Hi, be- 
have differently. For example, Fe ll does not respond in the same way to changes in CL2 and 
CL3. Figure m shows the great potential of the different IR lines to constrain the clumped 
structure of the stellar wind and demands the following detailed discussion. 

The general impact of clumping on line profiles that was described at the beginning 
of this section will occur provided the ionization equilibrium is on the "safe" side. We 
consider the "safe" region to be where the population of the next ionization stage clearly 
dominates over the one the line belongs to (i.e., Hll^Hl for the hydrogen lines). Noting, 
however, that ionization depends linearly on density whereas recombination is proportional 
to p^, a "changing" ionization situation may occur, where two adjacent ionization stages 
have similar populations. In such a case clumping, which enhances recombination, will cause 
a net reduction of the mean ionization. This will result in weaker lines. 

Finally, in the infrared, via bound-free and free-free processes oc p^), not only the 
lines but also the continuum will depend on clumping, resulting in high sensitivity of the 
continuum-rectified line profiles to CLi and CL2. 

One may, therefore, distinguish between lines formed on the "safe" region and those 
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arising from the "changing" region. Within the parameter domain of the two LBVs studied 
here we find that Hi, Sill, and Mgll hnes and also the Fell photospheric hnes are formed 
in "safe" regions, while He I and Fell] lines arise from "changing" regions. 

Increasing the clumping (decreasing the CLi value), or alternatively decreasing the 
velocity at which clumping sets in (decreasing CL2) results in stronger Sill and Mgll lines, as 
shown in Fig. |3]-left. The strong H I lines are formed further out than the Si ll and Mgll lines 
and their strengths should in principle show no sensitivity to CL2. However, the continuum 
is clearly affected by clumping. Thus the stronger the clumping (lower CL2), the stronger 
the continuum and the weaker the resulting line-to-continuum ratio, as clearly shown by B7 
and Ba in Fig. HJ-left. On the other hand, weaker Hi lines such as Hui4 or Bio and Bn form 
much closer to the photosphere and tend to brighten with increasing clumping. The weak 
Fe II lines formed close to the photosphere react in basically the same way as the continuum 
and thus the normalized spectra of them show no changes with changing clumping. Their 
near independency allows these lines to be used as Fe abundance indicators (see below). The 
Fe 11] lines, formed even beyond the H I lines, are affected by two competing processes. On 
one hand, increasing the extent of the clumped region (decreasing CL2) results in a reduction 
of the Felll/Fell ratio in the wind. Since Felll remains the dominant ionization stage in 
the Fell] line formation zone, this change will cause a slight increase in the strengths of 
the Fell] lines. On the other hand, as the continuum increases with increasing clumping, 
the line-to-continuum ratio decreases. Thus the two processes counter-balance (see Fig. H]- 
left). Finally, the He I lines, which form close to the photosphere, show weak continuum 
dependences, but high sensitivities to ionization/recombination. Thus, starting with the 
model with the highest CL2 values, the He I lines are not affected by clumping, but the 
stellar parameters produce strong ionization in the inner parts resulting in overly strong 
emission (He I 1.700 /im) and line filling (He I 2.112 /xm). However, as clumping is enhanced 
in the line formation zone, recombination starts to dominate over ionization and the He I 
line emission weakens, the lines are no longer as filled, and start to appear in absorption. 

Regarding clumping in the outer parts of the wind, it can be seen on the right side of 
Fig. m that only the strong Hi and Fell] lines react to CL3. Note that Ba which forms 
further out than B7 is more sensitive to clumping and the observed ratio of the two profiles 
may be used to determine CL3. For winds of significantly lower density, these lines will 
form further in and show little or no dependence on CL3 (e.g., FMM362). The Fell] lines 
are more sensitive to CL3, primarily due to the coupling of the Fe ionization structure with 
that of hydrogen thru charge-exchange reactions in the outer wind zones where H ll starts 
to recombine. Due to the difference in H and Fe abundances, a small and hardly noticeable 
change in the ionization of hydrogen will be amplified in the Felll/Fell ratio, resulting in a 
large change in iron recombination. Thus, decreasing CL3 dramatically enhances the Fell] 
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lines, as shown in on the right side of Fig. |H The H I hnes behave similarly to the way they 
do in "safe" regions. 

For the Pistol Star the best fits to the various line profiles suggest a decreasing clumping 
factor that becomes unity in the outer wind, whereas for the less dense wind of FMM362, the 
line profiles are best matched by a constant clumping factor. From the clumping estimates, 
we find uncertainties of 0.10 dex in the mass-loss rates. The value of CLi = 0.08 obtained 



for each st ar is somewhat low compared to the val ues derived 



(CLi=0.5,lNajarroJl2Q0ll) or AG Car (011=0.25, 



with those derived for WR stars (CLi~0.1, iHerald et al. 



Groh et 



al 



'or ot her LBVs such P Cygni 



2001 



20061) and is more con sistent 
Naiarro et aDl2004h . We 



note, however, that the LBVs have higher He abundances than P Cygni, pointing to a more 
evolved status, closer to the WR phase. 



4.3. H/He ratio. Breaking the degeneracy 



HiUier et all (119981 ) showed that for HDE 316285, an LBV-like star with similar Tes 
and slightly higher wind density than the Pistol Star, a degeneracy exists between the H/He 
ratio and the mass-loss rate. In principle, fits of virtually equal quality could be obtained 
with H/He ratios varying from 0.05 to 10 by scaling the mass-loss rate. Such a degeneracy, 
if present in the Quintuplet LBVs, would imply that if the H/He ratio falls below a certain 
value (H/He<2), the resulting metal abundances could be scaled down to obtain the same line 
strengths. Hence one could obtain only an upper limit on the metal abundances, breaking 
this degeneracy is crucial to understanding the evolutionary status of these objects. 

Because of the lower wind density of the Quintuplet LBVs and the sensitivities of some 
of the infrared lines to the stellar parameters, we are able to break the H/He degeneracy 
and obtain robust estimates of their He content. Due to the high degree of clumping found 
in both objects, the r = 2/3 radius, where T^fi is defined, is reached at considerably lower 
velocities than for classical LBVs. Thus, wind speeds roughly between half and one-third of 
the sound speed are found in the Pis tol Star and FMM362 while class ical LBVs have wind 
speeds well above the speed of sound (IHillier et al.lll998l : iNajarro Il200ll ). This enables quasi- 
photospheric absorption lines to form. The He I 2.112 /xm line is the key in breaking the 
degeneracy. This is shown for FMM362 in Fig. [5], which contains model spectra computed 
for H/He ratios ranging from H/He=5.0 to 0.75, with mass-loss rates and metal abundances 
scaled and the other stellar parameters fine-tuned to reproduce the observed profiles of other 
lines. The figure shows that while identical Hi and He I (He I 1.700 fim) line profiles (also 
for the rest of hydrogen and metal lines) are obtained for all H/He ratios considered, the 
absorption depths of the He I 2.112 nm and He I 2.15 /im lines react sensitively to the He 
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Fig. 5. — Breakdown of the H/He degeneracy in FMM362. Models with H/He ratios ranging 
from 5.0 to 0.75 provide identical Hi and strong Hel line profiles, but weaker Hel profiles 
can be used to determine the He abundance (see text) 
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abundance. Both lines show that the best H/He value must lie between 3.33 and 2.25, and 
we find a most likely value of 2.8 (see Table 1). Similar behavior was found for the Pistol 
Star, where we obtain H/He=1.5. 



4.4. Metal Abundances 



For th e purpose of discu s sing m etal abundances (see Table 2), we adopt the solar com- 
position of iGreyesse fc Noeld ( 199311. Alth ough their abu ndances have been recently revised 



(lAsplund et al 



2005 



AUende Prietd 120081) fbut see also IPinsonneault fc Delahayd (120061 )) 



they are the ones used by llglesias fc Rogers! (119961 ) to compute stellar interior opacities 
and adopted in the most recent evolutionary mo dels for massive stars with rotation from 
the Geneva group (IMevnet fc Maederl l2003l. 12005). and the P adova tracks used for cooler. 



less massive stars ( iGirardi et al 



200C 



lutionary tracks for massive stars (ISchaller et al. 



Salasnich et al. 



1992 



2000h. Previously p ublished evo 



Mevnet et al. 



tables calcul ated with solar composit ion from Anders fc GrevessX ( 1989 ). which differ signif- 



19941 ) used opacity 



i cantly from lGrevesse &: Noeld (119931 ) only in Fe (A(Fe/H)=7.67 vs 7.50 in lGrevesse fc Noels 
Jl993hlR and very sl i ghtly in the CNO ratios (A(C/H)=8.56, A(N/H)=8.05, A(0/H)=8.93 in 
Anders fc Grevessel (il989h vs A(C/H)=8.55, A(N/H)=7.97, A(0/H)=8.87 in lGrevesse fc Noels 
(119931 )). Nevertheless, we have also listed in Table 1, in parentheses, t he measured abun- 
dances with respect to the solar Fe values from I Anders &: Grevessel (119891 ). Si and Mg are the 
sa me in all evolut i onary models, and have been only slightly revised downward (~0.05 dex) 
by lAsplund et al.l (120051 ). Thus, the reader should note that current discussions found in the 
literature on the derived a-elements vs. Fe ratio may depend critically on the assumed Fe 
solar abundance. 

Iron. Two types of Fe ll lines are found in the spectra §. The first are the strong semi- 
forbidden lines, including z^F9/2-c'^F9/2 1.688 /im and z^F3/2-c'^F3/2 2.089 /im, that form in 
the outer wind and have small oscillator strengths (gf~10~^). The second are the weak 
permitted (gf~l) lines connecting higher lying levels, such as the 4de^G-5p®F lines near 
1.733 or 6p^D-6s^D at 2.109 /im, that form much closer to the photosphere. 

The permitted lines are more robust iron abundance indicators, having only weak depen- 
dences on other parameters, such us turbulent velocity. The strengths of the semi-forbidden 
lines depend on the accuracy of their weak gf values, the mass loss rate and the run of the 



2A(X/Y)=log[n(X)/n(Y)] + 12 

•^The forbidden [Fell] 1.677/xni line present in the Pistol Star is not included in our models 
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iron ionization structure in the outer wind, which is sensitive to the hydrogen ionization 
structure due to the strong couphng to the Fe/H charge-exchange reactions. Since a change 
in the run of the clumping factor in the outer wind regions modifies the ratio of recombi- 
nations/ionizations in hydrogen, the semi-forbidden hues are diagnostic of the behavior of 
clumping there. From Fig. [1] it can be seen that that our model is able to simultaneously 
reproduce both sets of lines, providing constraints on both clumping and abundancefl 

We obtain roughly solar iron abundances for both LBVs, with ±0.15dex as plausible 
uncertainties (see Fig. [7]). Our results are similar to A(Fe/H)=7.59 recently derived by 



Cunha et al.l (120071 ) from their analysis of a sample of luminous cool stars within 30 pc of 
the Galactic Center. Note in Table 1 that the Fe abundance ratio has significant uncertainty 
due to the uncertainty in the Fe abundance in the Sun. 

Magnesium. The strongest Mg ll lines observed in the H and K bands share the 5p^P 
level. Those lines with it as the upper level, the 2.13/14 /xm and 2.40/41 fim doublets (see 
Figs [Hand [2]) are much stronger than those with it as the lower level (H band lines), revealing 
that pumping through the resonance 3s^S-5p^P line must be a significant populator of the 
5p2p levels. Pumping through the 3s2S[l/2]-5p2p[3/2] 1025.968A transition is very efficient 
due to Ljj3 fluorescence. This was conflrmed in models in which we decoupled the 5p^P[3/2] 
and 5p^P[l/2] levels (see Fig. [6]), resulting in Mgll 2.13/14 fim ratios much higher than 
observed. 

The relevance of this process can be easily followed in Fig. El which displays the behavior 
of the doublet as a function of the choice of the Mgll atom and the turbulent velocity. The 
latter refers to the flxed Doppler width used in our models to compute the level populations. 
In the left panel of Fig. [6] the levels are considered to be decoupled (i.e., the number of 
superlevels in the model atom, NS, is set to the total number of levels in the full atom, NF). 
Because Ly/? lies closer to 3s2S[l/2]-5p2p[3/2] (Av=72kms-^) than to 3s2S[l/2]-5p2p[l/2] 
(Av=114kms~^), and because in these LBVS the terminal velocities and wind densities de- 
termine the line formation zones), only the 5p^P[3/2] is pumped thru fluorescence. Indeed, 
it can be seen that as the turbulence velocity is increased, the overlap between Ly/5 and the 
Mgll line increases, as does the population of the 5p^P[3/2] level and the strength of the 
Mgll 2.13 yum line increases, while the longer wavelength Mgll 2.14 /im line is unaffected. On 
the other hand, if the Mgll model atom has both levels combined into a superlevel (NSt^NF, 
Fig. [6]-right), the observed ratio is reproduced. Furthermore, increasing the turbulent veloc- 
ity, and hence the effect of fluorescence, increases the pumping of both levels equally and thus 



''We think that the shght missmatch m the Fell] 2.117 fim hnc in both objects is related to the accuracy 
of the f value. 
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increases the strength of the doublet with a constant ratio between its components. From 
Fig. [6] it can be seen that our assumed colhsion coefficients connecting the Mgll 5p^P[l/2] 
and 5p^P[3/2] levels may be too low. This comparison illustrates the importance of making 
the correct choice of model atoms for quantitative spectroscopic analysis. 

Due to fluorescence coupling, the Mgll K-band lines show a stronger dependence on 
turbulent velocity than do the H-Band lines. We estimate about twice solar Mg abundance 
and an associated uncertainty (see Fig. [7]) of about it0.25dex (due to uncertainties related 
to the fluorescence contribution). 

Silicon. The Sill doublet 5s^Si/2-5p^P3/2 1.691 fim and 5s^Si/2-5p^Pi/2 1.698 /im con- 
stitutes a powerful diagnostic tool, as it appears in emission for only a very narrow range of 
stellar temperatures and wind density structures, indicating the presence of amplified NLTE 
effects. However, since it forms at the base of the wind, its strong dependence on the details 
of the velocity field there hinders a precise silicon abundance determination. 

Instead, we use the well-behaved recombination line Sill 3s^6g^G-3s^5f^F at 1.718 fim 
which shows a stronger dependence on the silicon abundance. Once again, a realistic mapping 
of full- to super-levels in our model atom is required. From our model fits (see Figs. [1] and [2]) 
we derive roughly twice solar abundance (±0.20 dex) for silicon in each LBV, similar to 
magnesium (see Fig. [7]). 

Other elements. One might expect a number of oxygen lines might be detectable in 
infrared spectra of LBVs: i.e., strong Ol lines at 2.763 /im, 2.893 fj,m and 3.098 /xm and 
weaker lines at 1.8243 fim, 3.661 /im and 3.946 /zm. Several of these, but not all, are prob- 
lematical from ground-based observatories. Unfortunately our data set only encompasses 
the Ol 1.745 /im line which is blended with a stronger Mgll line. Thus, we defer an at- 
tempt to estimate the oxygen abundance until high resolution observations of unblended 
lines can be obtained. Determining the oxygen abundances in these objects will provide 
crucial constraints on their evolutionary status. Models show that when H/He<1.50 oxygen 
has reached its maximum depletion within CNO equilibrium, while a significantly higher O 
content should be present on the stellar surface for H/He values around 3. The results in 
Table 1 then predict that the Pistol Star and FMM362 have different oxygen abundances. 
On the other hand, if an LBV has a H/He<1.50 but is still hydrogen rich, the oxygen abun- 
dance determination, expected to be ~0.04 of the original value, will provide a measure of 
the metallicity of the natal cloud. High resolution L-Band spectra of the Pistol Star should 
be able to address this issue. 

The only detected sodium lines are the well-known doublet at 2.206/9 /im, from which 
we obtain a very high abundance, ~ 20 x solar. The strong observed emission of this doublet 
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in the K-band spectra of other LBVs has been noted previously by iHiUier et aLl (119981 ) . 
Interestingly, our models display only a minor dependence of these lines on clumping. On 
the other hand, the strengths of the sodium lines might not indicate extraordinary sodium 
abundance if the lines are produced by fluorescence of circumstellar material, a component 
that we do not model. 



5. Discussion. 



Our results suggest solar Fe abundances and approximately twice-solar a-element abun- 
dances for the Quintuplet LBVs. Presumably, these abundances were the same in the gas 
that condensed to form these stars and the other stars in the Quintuplet cluster and indeed 
in the whole of the present-day Galactic center. The results can be discussed in the context 
of similar measurements of Galactic center objects and with respect to the trend one might 
expect if the region is an inward extension of the disk or the bulge. In addition, the ratio of 
Fe to a-elements might be used to decipher the star formation history in the Galactic center. 

Table 2 displays a number of recent determinations of stellar metal abundances in the 
Galactic Center together with the above mentioned three reference patterns for solar abun- 



danc e s. The values derived for Fe abundances in cool stars a gree with our result (ICarr et al. 



2000l : iRamirez et al.l 119971 . Il999l . l2000f ). ICunha et all (120071 ) find a very narrow range of Fe 
abundances clustered around the solar value f or a population of cool stars in the central 
30 pc. Of particular interest is star VR5-7 from lCunha et al.l (120071 ) sample which is located 
in the Quintuplet Cluster and shows A(Fe/H)=7.60 and A(Ca/H)=6.41 

There are relativ ely few measurements of the a-element abundances ([a/Fe]) in GC stars. 
Najarro et al.l (120041 ) find solar abundances (as defined in this paper) for hot stars in the 
Arches cluster based on the oxygen abundance and, to a lesser degree, carbon ab undance, and 
adopting the canonical solar value of A(0/H)=8.93 (I Anders fc Grevesselll989l ). Those esti- 
mates assume that nitrogen has reached its maximum surface abundance value. Evolutionary 
models in dicate that 95% of t hat value is already attained by the time that H/He<2 (by 
number). iNajarro et al.l (120041 ) followed the metallicity patterns from the Geneva evolution- 
ary models and assumed no selective enrichment of CNO or a-elements vs Fe, in concluding 
that the stars in the Arches Cluster have solar a-element abundances. However, estimates of 



solar abundances have varied considerably over the past 15 years (e.g. lAUende Prietol 12008 



, see also Table 2). Th us, depending on the assumed solar CN O composition, the derive d 
nitrogen ab undance bv iNaiarro et al.l (120041 ) could imply solar (I Anders &: Grevessd Il989l ) , 
1.2 X solar (IGrevesse fc Noelslll993l ) or 2.0 x solar (lAsplund et al.ll2005l ) CNO composition. 
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Table 2. Metallicity Estimates in the GC 



A(Fc/H) 



A(0) 



ECNO 



Z(N) 



A{Mg/H) A{Si/H) A{Ca/H) 



Sun 

Anders c 



Grevesse (19891 7.67 
7.50 



Grevesse & Noels fl993) 
Asplund ct al. (2005) 



7.45 



;.93 

;.87 
;.66 



1.6 
1.4 
0.82 



1.5 
1.3 



7.58 
7.58 
7.53 



7.55 
7.55 
7.51 



6.36 
6.36 
6.31 



Hot Stars 
This work 
^Jaj^ro_et_aL 
^eballe_et_aL 

Martins et al. 
Martins et al. 




7.54±0.15 



;.91±0.12 



1.57±0.45 



1.43±0.42 
1.80±0.50 



7.85 



7.84 



Cool Stars 
Carr et al. (20001 
Ramirez et al. (2000) 
Cunha et al. (2007) 



7.48±0.13 
7.61±0.22 

7.59±0.06 9.04±0.19 



1.14±0.55 



6.71±0.14 



Note. — Compilation of estimated stellar metal abundances at the GC. All abundances are given as A(X/H) except column 4 
which displays the total CNO mass percentage and column 5 displaying the nitrogen surface abundance percentage Z(N). For 
the first two rows (solar values), column 5 gives the maximum nitrogen surface mass abundance, Z(N)max, reached during 
evolution according to the Geneva models. For the other rows, the average of the derived Z( N) for the s ample in each work 
is displa yed. All studies of h ot st ars make use of CMF GEN, while those of cool stars utilize IPlej lll992|) models grid. The 
resu lts oflNaiarro et al.l ll2004h and lMartins et al.] ll2007bl) are for stars in the A rches Cluster while tho se of lcteballe et al.l ll2006h 
and lMartins et al. I l l2007al ~ refer to the Central Parsec Cluster. Star AFNW in iMartnnset aL I ll2007al) has not being considered 
due to the poorer S/ N and resolution quoted by the authors. Th e results of Carr et al. ril2000h are for IRS7 while those of 
iRamirez et al.l 1 I2OOC1I ) are for six M supergiants and three giants. ICunha et al.l ||2007| ) extends the sample of iRamirez et al] 
hood) and also computes CNO and a-elements abundances. 
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Recently iMartins et al. J2007al lbh have analyzed a larger sample of hot stars in the 
Arches and Central Parsec clusters and find similar results (see Table 2). Interestingly, if 



one considers only the objects in iMartins et al.l ( 2007bl ) with He/H> 0.1 and those with 
Z(C)< 0.05, i.e., fulfilling the condition to be close enough to Z(N)max, the average value of 
Z(N) is 1.7. 



Geballe et al.l (120061 ) estimate roughl y solar oxygen abun dance, A(0/H)=8.91, in IRS 8 



an Olf supergiant near the central parsec. ICunha et al.l (120071 ) find <A(O/H)>=9.04 ([0.37]) 
and <A(Fe/H)>=7.59 ([0.14]) for their sample of cool stars, where the numbers in brackets 
are the ratio with respect to the solar value ii i dex. This i r nplies [O/Fe]=0.22, i.e. a clear 
enhancement over the solar ratio. Again, the ICunha et al.l (120071 ) measurements could be 
interpreted as indicating solar ratios in O over Fe if the solar O abundance in evolutionary 
models is used. It is crucial to have accurate solar abundances, and that values used in stellar 
evolution calculations should be consistent with these. An excellent example is attempting 
to determine whether the possible oxygen enhancement is due to a top-heavy IMF favoring 
ct-elements vs Fe enrichment, or simply an overall CN O and metal enhan cement. Thus, 
taking the CNO abundances for the GC objects from ICunha et al.l (120071 ) and assuming 
C/N equilibrium values one can interprete their results either as solar CNO with mildly 
enhanced (30%) oxygen (lAnders &: GrevessdllQSQl ) or a clearl y supersolar eiiyiron ment with 
a factor of 1.7 enhancement for C and N and 2.5 for oxygen (lAsplund et al.ll2005l ). 



Fortunately, there are other a-elements whose adopted solar abundances have suffered 
basically no major revision. Thus, we believe that the enhanced values obtained in this work 
for Mg and S i , roug hly a factor of two solar, together with the enhancement of Ca found by 
Cunha et al.l (120071 ). are a strong indication of the enrichment of a-elements compared to 
Fe. 



Our results run counter to the trend in the disk ( Rolleston et. al.ll2000l : ISmartt et. al. 



Martm-Hernandez et al. 


2003), and are 


more 


(Froeel et al. 


1999 




Feltzine & Gilmore 


2000) 



disk does not extend inward to the GC, so that material is dragged into the central molecular 
zone from the bulge rather than from the disk. Another possibility is that the GC stars are 
forming out of an ISM that has an enrichment history distinctly different from that of the 
disk. At this point, further studies of the a-elements vs Fe would be useful. Future high 
S/N and high resolution spectroscopy of the Ol lines in LBVs and K-band spectroscopy 
of WNL stars in the same cluster (Najarro et al. in prep.) will provide two independent 
measurements of the original oxygen content, and thus set definite constraints on metallicity. 

The modest enrichment in a-elements versus Fe that we find in the two Quintuplet 
LBVs is consistent with a top-heavy IMF in the GC (iFiger et al.lll999al ). In such a scenario. 



enhanced yields of a-elements compared to Fe are expected through a higher than average 
ratio of the number of SNII vs SNIa events (IWheeler et al.lll989t ICunha et al.l 120071 ) . 
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Fig. 6. — Influence of Ly/? fluorescence and choice of Mgll model atom on tlie Mgll K-band 
lines. 



-27- 




X 



tS] 



as 0.9 : 

6 

f-, 
o 
2 



0.8 : 



1.2 : 

1.1 - 



1 .0 "'^-^ 



Mgll Hel 



Mgll Br^^ 




1.50 
1.40 
1.30 
1.20 
1.10 
1.00 



1 

Mgll 


1 

Mgll -. 




Fell 


:| 






i 1 














1 1 



1.673 1.676 1.679 1.682 2.1350 2.1393 2.1437 2.1480 




1.10 - 



1.05 - 



1.00 ' 



0.95 - 



1.7160 1.7177 1.7193 1.7210 1.6960 1.6970 1.6980 1.6990 

Wavelength 

Fig. 7. — Error estimates of Fe (upper-panel), Mg (middle-panel) and Si (lower-panel) abun- 
dandes. Dashed lines (red) correspond to our best model fitting the observed (black-solid) 
diagnostic lines of FMM362. Long-dashed (green) and dashed-dotted (blue) lines correspond 




